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0\ , Abstract 

!>; 

■r^ij- I We consider the linear perturbations for the single scalar field infiation model interacting 

with an additional triad of scalar fields. The background solutions of the three additional 
scalar fields depend on spatial coordinates with a constant gradient a and the ensuing evolu- 
tion preserves the homogeneity of the cosmological principle. After we discuss the properties 
p\ , of background evolution including an exact solution for the exponential-type potential, we 

C^ . investigate the linear perturbations of the scalar and tensor modes in the background of the 

slow-roll inflation. In our model with small a, the comoving wavenumber has a lower bound 
~ aMp to have well-defined initial quantum states. We find that cosmological quantities, for 
instance, the power spectrums and spectral indices of the comoving curvature and isocurva- 
ture perturbations, and the running of the spectral indices have small corrections depending 
on the lower bound. Similar behaviors happen for the tensor perturbation with the same 
lower bound. 
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1 Introduction 

Inflation [Tl|2] has been so far very successful to account for the observational data [3lll]. It is based 
on the simple idea that the Universe went through a huge accelerating expansion during the early 
stage of its evolution driven by a single (multiple) scalar fleld(s). It provides not only spatially 
flat and homogeneous universe but also a source for small primordial perturbations which are the 
origin of the large-scale structure of our Universe today. Even though the data are compatible 
with the single fleld inflation, there have been considerable interests in multi-fleld inflation. 

From the point of view of cosmological perturbations [5] , multi-fleld inflation exhibits distinc- 
tive features of the non-Gaussianity ^Ul and generates isocurvature (entropy) perturbations in 
addition to adiabatic perturbations of the single scalar fleld model. There exists an extensive 
amount of literature how these aspects are incorporated into the speciflc multi-fleld inflationary 
models p-[TU]. In particular, the presence of isocurvature modes corresponds to relative pertur- 
bations among the various matter flelds and affects the flnal curvature perturbation at the end of 
inflation by acting as a source term in the evolution equation for the curvature perturbation [TO] . 

So far, most of analysis focused on the time- dependent backgrounds in which each of the 
scalar fleld is a function of time only. However, growing interest in this fleld motivates to seek 



other theoretical possibihties. One of the alternative approaches is to consider spatially dependent 
backgrounds [H]. It seems that relatively little attention has been paid to spatially dependent 
backgrounds in inflation. This might be due to the obstacle which is inherent in these solutions 
in cosmology: in general, these configurations are not compatible with the cosmological principle 
of homogeneity and isotropy. However, the confiict can be avoided in the nonlinear sigma models 
which include a triad of scalar fields, a" (a = 1, 2, 3). It is known that by correlating each scalar 
field with spatial coordinate, a" ~ x'^, one can maintain the isotropy and homogeneity in a wide 
class of field theory models whose Lagrangians are functions of X = —\g^^ hab{cr)d^a°'dya'^ without 
the potential for the a"" fields. 

The ansatz for scalar fields of a nonlinear sigma model, o"™" ~ x™" {i = 1, ■ ■ ■ , A^) with the 
extra dimension A^, first appeared in the higher dimensional gravity theory and the solutions 
of the scalar fields trigger spontaneous compactification of the extra dimensions [T21IT3] . The 
compactified extra space is isomorphic to the target space of scalar fields. A similar ansatz which 
break the diffeomorphism invariance was used to give masses to gravitons as a Higgs mechanism of 
gravity [13] • Therefore, the coordinate dependent ansatz in the nonlinear sigma model in gravity 
theories can be considered as a method of constructing the massive gravity theories [15]. Recently, 
the four scalar fields were combined into de Sitter target space and used in describing the late-time 
accelerating universe [16]. See also [TJIIIH] for related topics. 

Motivated by the nonlinear sigma models inherited from the higher dimensional gravity theo- 
ries, in this paper, we consider a nonlinear sigma model with a potential as an action of matter 
fields in four-dimensional gravity theory, 

S = j d'^x^^g (^R-^g'^'^h^,^X^)d^'^^d,a^ - V{a)\ , (1.1) 

where m, n = 1,2,3,4 and Mp denotes the Planck mass, Mp = (8nG)~^^'^. Here hmn can be 
considered as an internal metric of a four-dimensional Riemannian manifold. To construct an 
infiationary model, we consider a specific choice for the internal metric and potential, 

Kb = f{v)Sab, haA = 0, /l44 = 1, 

Via) = Vi^), (1.2) 

where (^ = a^ and f{f) is an arbitrary positive function of yqj. Then the resulting model can be 
considered as the single scalar field model interacting with a triad of scalar fields a", where the 
coupling among the infiaton (f and the triad of scalar fields is non-minimal, ~ f{(f)X. 



^It is to be pointed out that the four fields have all positive kinetic energy unlike the previously mentioned cases 
of massive graviton and de Sitter target space where the scalar field if has a negative kinetic energy. 



We try to describe the cosmological inflation in terms of tlie action (11. ip . As a background 
solution with the choice of (II ■2p . we apply the ansatz a"" ~ x" for the background evolution, 
assuming that other dynamical fields depend on the cosmic time only. As we discussed previously, 
the resulting cosmological evolution becomes homogeneous and isotropic with the FRW metric. 
For the background evolution, we find an exact solution describing the power-law inflation and 
try to figure out the properties of the slow-roll inflation of our model. 

The main purpose of this paper is to investigate the cosmological perturbations in this modelo 
We consider the linearized scalar and tensor perturbations for the slow-roll inflation with a small 
spatial constant gradient a for the triad of scalar fields|j The scalar modes have two degrees 
of freedom, one from ip and the other from the scalar mode of the triad. They are decoupled 
from each other in the minimal coupling (/(v?) = 1) and have the form of Sasaki-Mukhanov 
equation [12] in the spatially flat gauge. We obtain the leading contributions of a for the power 
spectrums and spectral indices of the comoving curvature and isocurvature perturbations. 

One interesting feature of our approach is that in the slow-roll case, the requirement of unitarity 
for the initial quantum state imposes a lower bound on the comoving wave number. The existence 
of this limit is a consequence of the nonvanishing aforementioned free parameter a associated with 
spatial condensations. 

The paper is organized as follows. In section 2, we start from the action (II. ip with the choice 
(II. 2p . Using the ansatz a"" ~ x"" for the triad of scalar fields, we obtain an exact solution for the 
power-law inflation and investigate the properties of the background evolution for the slow-roll 
inflation. In section 3, we write the perturbed equations for the scalar and tensor modes and 
discuss the gauge conditions for the scalar perturbation. In section 4, we describe the behaviors 
of the perturbations up to leading order of a and obtain cosmological quantities. We conclude in 
section 5. 

2 Background Evolutions 

We start from the action (II. ip with the choice (II. 2p . 



S = d'^Xy/^ 



M|„ 1 „„„ „ 1 



R-- g^'d.ifd^cp - -f\v)9"'d,a'^d,a^ - V(^) 



2 2 '^^ ^ 2 



(2.3) 



^ Cosmological perturbations with the triad of scalar fields having solutions a"^ ^ x"" were considered before 
in [11| . The main interest of the paper was to create statistically anisotropic and inhomogeneous perturbations 
to account for certain anomalies in the temperature anisotropics in the CMB. The metric perturbations were not 
dealt with in pTj . 

■^The three perturbed modes of the triad of scalar fields are decomposed into one scalar and two vector modes. 



where a^'s have an S0(3) symmetry. This action can be considered as the single field infia- 
tion model with a triad of scalar fields having noncanonical kinetic terms. We read the energy- 
momentum tensor as 



T^y = d^ipd^ip + f{ip)d^a°'d^a'' - g^^ 



\f\v)9''^d^a'^dpa'^ + ^ g'^^d^ipdp^ + V{^) 



and equations of motion of the scalar fields ip and a"" as 



(2.4) 



(2.5) 



where V^p = dV/dip. 

The background FRW metric is given by 

ds^ = -dt^ + a{tf{dx^ + dy'^ + dz^) , (2.6) 

where a{t) is the scale factor. In this paper, we consider an ansatz for the scalar field a" as 



a" = M^a x" 



(2.7) 



where a is an arbitrary dimensionless constant. For the similar spatially dependent ansatz for 
S0(3) symmetric fields, see also [H]. Under the assumption that the scalar field ip is spatially 
homogeneous, we easily see that the ansatz (12. 7p satisfies the equation of motion for a'^ in (12. 5p . 
Remaining equations of motion for g^^ and ip are given by 

P 



H = 



3M2' 
P + P 



2M2 ' 



(2i 



3M^a 



4^,2 



ip + 3Hip+^^f^ + V^ = 0. 



where f^ = df/dip, H = d/a, and 



1 , 3M^a2 
^ 2^ 2a2 ■' 



(2.9) 



As we see from these background equations, the homogeneity of the background evolution is not 
spoiled though we take the spatially dependent ansatz for a"- in (12. 7p . The /-dependent terms in 



(12. 9 p are originated from the contributions of a"". Therefore, the equation of state of a- fields is 
Wa = —1/3. For concreteness, we fix the function f{ip) as 

/(^) = e2«^/^^p, (2.10) 

where ^ is an arbitrary constant. In the case ^ = 0, i.e., /(v?) = 1, the role of the a^-dependent 
term in (12.81) is the same as the curvature constant determining the spatial curvature in the 
Friedmann equations. Since we are considering the positive a^, the background evolution under 
the spatial dependent solution (12. 7p corresponds to that of the open universe. However, our model 
is different from the single field inflation model on the background metric of the open universe. 
As we will see later, there are nontrivial roles of a" in the perturbation level. 

The main purpose of this work is to investigate the behaviors of our model (12.30 from the point 
of view of linear perturbations. Before we move on to the subject of the perturbation, we describe 
the background evolutions of our model by considering two inflationary scenarios, the power-law 
inflation and the slow-roll inflation. 

2.1 Power-law inflation 

It is well-known that under an exponential-type potential the single scalar field model has an 
exact solution describing the power-law infiation [20]. Comparing with the single scalar field 
model, there are a-dependent terms in (12. 8p . which are originated from the spatially dependent 
background solution for a"'. As we will see in this subsection, the exact solution for the power-law 
infiation in the single scalar field model can be a solution of the equations of motion in (12. Sp with 
some deformations by parameters, C, and a. 

As in the case of the single scalar field model, we consider an exponential-type potential 

y(v9) = K)e-^^/*^p, (2.11) 

where Vq and A are arbitrary constants. We consider an ansatz which describes the power-law 
infiation, 

a(t) = ao(Mpt)'^, ^{t)=Mp(hn{Mpt)+ifo] , (2.12) 

where oq, n, and (po are dimensionless constants. 

When a = 0, the background equations reduce to those of the single scalar field model. The 
constant parameters in the power-law solution ( ]2.12p are determined as 

"=^' ^° = ^" (2M^|^) ' ^2.13) 



and ao becomes a free parameter. 

However, in the case of a 7^ 0, all constant parameters of the power-law solution (12.12^ are 
fixed by 



n = l + -, 



«2A2 / AVn ^'^/^ 



"° 2(A2 + 2Ae-2) V2M^(6e + 3Ae + 2)J ' ^^'^^^ 

c/^o = - In 



A V2M|(6e2 + 3A^ + 2), 

The positive definiteness of the right-hand side of the second condition in f l2.14p restricts the value 
of A, 



A > -e + Ve' + 2 or X<-^-^/eT2. (2.15) 

In order to have an accelerating unverse, we also have a constraint, f > 0, from the first relation 
in fl2.14p . Therefore, there are two possibilities for the parameters, A and C,, i.e., A > & .^ > 
or A < & .^ < 0. In the case A>0&.^>Owe have the first relation in f l2.15p . while in the 
A < & ^ < case we have the second one in fl2.15p . Then the range of n for both cases can be 
expressed in terms of ^ as follows, 

i< n < 1 + e' + ve^ + 2^2 . (2.16) 

In the case ^ = 0, n is fixed to unity. Therefore, there is no inflation in this case. In this paper, 
we consider A > & ,^ > case. 

The power-law solutions in 02.13P and 02.14p are special solutions of the background equations 
of (12. 8 p since we started from a special ansatz in obtaining the solutions. For this reason, it is 
understandable that one cannot obtain the solution (I2.13P by taking a — )■ in (I2.14p . That is, 
these two solutions satisfy different initial conditions in the a — )■ limit. The general solution 
under the exponential-type potential (12. lip for the single scalar field model was obtained in 



2.2 Slow-roll inflation 

Next, we consider the slow- roll inflation with a potential satisfying the slow-roll approximations, 

02/2 < ^, <^ < SHcp. (2.17) 

To reflect these approximations, we introduce the slow-roll parameters, 

V^^,- (2.18) 



0^ ^V, 



ipj 




0.02 0.03 

a 



Figure 1: Plots of e-foldings for the quadratic potential V = ^m ip on the parameter space of a 
and ^. We choose the initial condition as qq = 1, V^o = 16Mp, and 0o = with m = 10~^Mp. 



With the approximations in (12.171) . the background equations of motion of (12. 8p reduce to 

1 



~ 2a2 



/ 



3M2"^' 



3H^ + ^^U + V,^0, 



(2.19) 



where f^^ = jj-e^^'^'^^^ for the case (12.101) . The background dynamics of ^p is governed by the two 
papameters a and ^ as well as the shape of potential V{ip). 

We can compute the number of e-foldings using the slow- roll equations of motion of (I2.19p . 



M 



Hdt 



ClLp 



Mi 



,. y a'M^f 



+ 



2a2 



P^.e }^+f^!ii^ 



dip, 



(2.20) 



where (pi and (pe are the values of ip at the initial time and at the end of inflation, respectively. 
For a fixed value of a, the number of e-foldings decreases with increasing ^ since the .^-dependent 
term only appears in the denominator in (I2.20p . For concreteness, we consider a massive scalar 
field with a potential 



V{(p) = -mV^- 



(2.21) 



Then the number of e-foldings in fl2.20p is given by 



For a small value of a, we can expand the number of e-foldings as 

^. 1 r,[l-(^-^)^ + 0(a.*)]#. (2,23) 

In this limit, the a-dependent term in fl2.23p is also small and the contribution of this term to the 
number of e-foldings is negligible and independent of the magnitude of $,. However, the sign of 
^ — Mp/(f determines the negative or positive contributions to the number of e-foldings. Therefore, 
there is a critical point for the value of ^ near 

^ ^ Mp/(f. (2.24) 

The scalar field (f is varying as the cosmic time evolves. However, it is almost a constant at the 
early stage of the inflation, according to our assumption of the slow-roll inflation. For this reason, 
the value of ^ in fl2.24p can be the critical point. The behaviors of the e-foldings for the quadratic 
potential fl2.2ip are given graphically in Fig.l. 



3 Linear Perturbations 

In this section, we discuss the linear perturbations for our model. We introduce a new gauge 
invariant variable in the scalar perturbation, which is originated from the space dependent back- 
ground solutions of the triad of scalar fields, a". After we take the spatially flat gauge which is 
convenient in our later analysis for the linear perturbation, we write the general form of the comov- 
ing curvature perturbation and the isocurvature perturbation. We also discuss the perturbation 
for the tensor mode. 

3.1 Scalar mode 

In this subsection we consider the linear scalar perturbation for our model (12. 3p . The linear scalar 
perturbation of the metric is given by 

ds'^ = - (1 + 2A) dt^ + 2a diB dtdx' + a^ [ (1 - 2^^) 5ij + 2didjE\ dx'dx\ (3.25) 



where A, B, ip, and E are four scalar perturbation modes. We also consider perturbations of the 
scalar fields, 

V9(t, x) = (p{t) + 6(p{t, x) , 
a'^(t,x) =a"(x) + 5a"(t,x). (3.26) 

Before we write the perturbation equations for scalar modes, we discuss the contribution of 
(5(T* to scalar modes. At first, we read the (O,i)-component of the perturbed Einstein equation 
5G^i = Mp 2 5T0 as 

-2 di (ha + tp] = Mp^ST^i = Mp 2 (-diiipSifi) - afSa' + ^diBj . (3.27) 

Taking the curl of both sides of (I3.27p . we obtain 

e^jkd^Sa^ = 0. (3.28) 

From this relation, we see that Sa"^ appearing in the linear scalar perturbation equations has no 
contribution to the perpendicular mode. For this reason, we can seel 

Sa' = ydiU, (3.29) 

k 

where k is the comoving wave number and u is the perturbed scalar mode. Here we introduce 
1/k factor in order to define the canonical kinetic term for the scalar mode u in the perturbed 
Lagrangian. Then we can write the variation of the (0, z)-component of the energy momentum 
tensor as 

5r° = diSq, (3.30) 

where the scalar part of the three momentum Sq is given by 

Sq = -^S^ - -^fu + ^—fB. (3.31) 

k a 



*In general, Sa^ is decomposed as Sa"^ = jdiU + 5a]_ with the perpendicular mode Sa]_ satisfying di5a]_ = 
in the linear perturbation. When we consider the perturbation of the vector mode, Sa\^ can contribute to the 
perturbed equations. However, we do not consider the perturbation of the vector mode in this paper. 
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3.1.1 Perturbation equations 



Inserting the relations for perturbed modes, fl3.25p . (13.261) . and (13.291) into the Einstein equation, 
we obtain the linearized perturbed Einstein equations SC^^^ = -^ST^^^: 



{0,0) -.SHi HA + ij)--V 



^ + H{a'^E-aB 



Sp 



(0,0 -.HA + ij 



5q 

2M2 



iIj + 3H(HA + '^) +HA + 2HA + -^V^D 



5p5V + nv 



2M2 



5v - -,dAD 



Ml 



(3.32) 



where 



dV M^a^ I „/ n 



E]+'iiP + 






Mta 



4^,2 



d 



u 



3a2 ^ V ^akMl -^J+'^+ mTJ 



D = A-i)-H[a^E ~aB] ~—[a^E-aB] , 



W 



2Mt 



a 



f 



^^^^■(^ 



-E] --V 



u 



3 ■ \akMl 



-E 



(3.33) 



From (12. 4p and (12. 5|) . we obtain the perturbed equations for 6ip and u without gauge fixing. 



1 



2M'^a'^ / V^ 



50 + 3H6^ - -V'6v + "^'^i; V (^-^ + 3^5^ ) + ^6^ 



a^ 



2A(p + <oEA^ + 



.2/i? 



A + 3V^ + V" E 

V a 



2M|a2 _c 



/ 3^-V2i? 



i2+|3i/+|^,n 



1 „2 «^^P 
_V\ = ^ 



,B + 2,|^ + || B + ?|^H.(^_^._V^S) 



(3.34) 



Using the relations in (I3.3ip and (13.33^ . we can write the comoving curvature perturbation j22t 
and the isocurvature perturbation for multiple scalar fields [TU] as 



TZ = ^ip 



H 



P + P 



P P 



(3.35) 
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3.1.2 Gauge condition 

As we see in f l3.29p . we have one additional scalar degree of freedom due to the contribution of 
5(t'^. Therefore, we have six degrees of freedom for the scalar perturbation modes. Eliminating 
the gauge degrees of freedom and imposing the constraints of the Einstein equation, we have two 
physical degrees of freedom. Since we have considered the spatially dependent background solution 
for cr'^, we have slightly different gauge invariant quantities from the well-known multi-field models. 
Under the change of coordinate, 

x^ ^ x^ + /3^, (/3° = /3, /3^ = 5*^9j7), (3.36) 

the scalar modes of the perturbed metric (I3.25p transform as 

A -^ A-/3, 

B -^ B + --a^, 
a 

E -^ E--f. (3.37) 

The scalar modes of matter fields also transform as 

5(p — )■ 5(p — [5^d^ip = 5(p — (pl3, 
da" -^ 5a^ - P^'d^a" = jS^'^dbU - aM^S'^'^db^, (3.38) 

fv 

where we used the background solution a" = aMpx"' of (12.71) in the last step. From (I3.38P we can 
read off the gauge transformation of the scalar mode u as 

u ^ u- kaM^j. (3.39) 

From these transformation rules, we can obtain several gauge invariant quantities. For later 
convenience we define those quantities, 

d 



^ = ij + aH{B-aE), 



Q^ = 6(p- —ijj, 



H' 

Qu = u- akMpE. (3.40) 



The new quantity Qu was introduced due to the spatial dependence of the background field a"'. 
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To solve the perturbed equations (13.321) and (13.261) . we take the spatially flat gauge {if) 
& -E = 0). In this gauge the perturbed equations are reduced to 






Q^ + 3HQ^+ — 



. +v,..5M^/ + H^^/IQ 



a^ 



am 



v^ 






g.+ |3i. + |^|g.+ ^^ 



k' 2a'Mlf\ 



2^akMp/^ ^ \ 

-^-^(Q^ + aBipj =0, 



(3.41) 



where Q^ = 6ip, Qu = u. The scalar modes and A and B satisfy the constraints which are given 
from the first two equations of (I3.32p . 



?,AH^ 



2AH 






0Q^, (a^ _a^B 



a 



I/. 



(3.42) 



The comoving curvature and isocurvature perturbations of (I3.35P are written in the spatially 
flat gauge as 



n = H 



S = H 



^Q,-aMlf{^^-<^) 



^ 



a? J 



ifiQ^ - Aif) - V^Q^ + -^f 



"^^P f I kQu 3?Q, 
T 



ifiQ^ - Aif) + V^Q^ - ^^/^ ^]vS A/p 



T 



^(»3-vy + 2!^/(j/-|j) 



«i(#H-V.)-2^/(ff-ft) 

(3.43) 



In the case a = 0, the quantities TZ and S are reduced to those of the single scalar field [TO] . 



3.2 Tensor mode 

In this subsection, we investigate the linear perturbation of the tensor mode for our model on the 
background described by the equations of (12.71) and (12.81) . Since the tensor mode is decoupled from 
the scalar and vector ones in the linear perturbation theory, we consider the perturbed metric in 
the conformal time to describe the tensor perturbation. 



ds"^ = a'^{T) [-dr'^ + {5ij + hij)dx'dx^) 



(3.44) 
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where dt = adr. The tensor modes hij satisfy the following conditions, 

dihij = 0, ha = 0. (3.45) 

Since the tensor modes h^j have two degrees of freedom, which are identified with the two polar- 
izations of gravitational waves, we expand the tensor modes as 

Mr.f)^-|jE^4. (3.46) 

where e^- is the polarization tensor satisfying the orthogonality condition e^e^"*, = 6y. Inserting 
fl3.44p and (13.461) into the action (12.31) on our background and solving the equation of the motion 
of fix, we see that the Fourier mode of fix satisfies 

fiUr-k)+(^e-^^fix{T-k) = 0, (3.47) 

where ' = d/dr and fix satisfies the normalization condition, 

/^A/i'I - ^>'a = «• (3.48) 

Using the relations for the background evolution in (12. Sp . we obtain a general differential equation 
for the tensor mode for our case, 

/^A + {k' - \MWf + ^(^'^ - 4aV))/iA = 0. (3.49) 

4 Perturbations in the Slow Roll Inflation 

As we see from the perturbed equations (13.411) and (I3.42p in the spatially flat gauge, the perturbed 
modes Q^p and Qu are completely decoupled in the a — )■ limit. Moreover, the resulting curvature 
perturbations, TZ and S in (I3.43p . are independent of Qu- Therefore, in the a — )■ limit, the linear 
perturbation of our model is the same as that of the single scalar field model. For this reason, in 
this paper, we only consider the case with nonvanishing a. 

When we consider the linear perturbations with nonvanishing ^ and a, in general the modes 
Q^p and Qu are not decoupled. However, in this coupled case it is very difficult to obtain some 
cosmological quantities, such as power spectrums and scalar/tensor spectral indices. We leave this 
general case for further investigation. In this paper, we restrict the linear perturbations to the 
case of the slow-roll infiation at a small value of a with vanishing ^. As we will see in this case, 
the modes Q^, and Qu are decoupled. However, TZ and S still have dependence on both Q^ and 
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In this paper, we will not consider the perturbation for the power-law inflation. The reason 
is the following. For the case of a 7^ with ^ = 0, the background equations in (12.80 have a 
special solution given in f l2.12p and f l2.14p for the exponential- type potential (12. lip . In this case, 
the solutions are reduced to 

^ X V 2Mp y ' ^ ^ V2A2 - 4 ^ ' 

The scale factor is a linear function of the cosmic time t, and so accelerating universe is not 
allowed. Presumably, if we can obtain the general solution like the work for the single scalar with 
exponential-type potential in [21], we may have an inflation era even in the case ^ = 0. However, 
it is beyond the range of this paper. Due to this reason, we do not consider the cosmological 
perturbation for the power-law inflation. In this section, we concentrate on the perturbation for 
the slow- roll inflation described in section 2. 

4.1 Scalar mode 

In the slow-roll approximation at a small value of a with .^ = 0, the perturbed equation fl3.4ip are 
reduced to 




Q^ + 3HQ^+{- + ^-^ + V^^\Q^c^0, 



Qu + 3i7Q„+ - + ^ Q„ - 0, (4.51) 

\ Ci (Jj I 

where we keep the leading order in slow-roll parameters of (I2.18p . Changing the cosmic time t 
into the conformal time r and introducing the Sasaki-Mukhanov variables, 

V = aQ^, U = aQ^, (4.52) 

we obtain the Sasaki-Mukhanov equations from (I4.5ip . 



K+[ kl-^^\U, = Q. (4.53) 



r2 
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The parameters in fl4.53p are given by 



.2_,2 «'M2 



rCi — /b 



6 



.2 _ , 2 , ll^^^g 

ft, 2 — ft- "I- „ ; 

/^i - 2 + 3e - ?7, 
3 

/X2-- + 6. 



(4.54) 



The /ci and k2 can be considered as the effective wavevectors for the Vk and hlk modes, respectively. 
In driving the Sasaki- Mukhanov equations in (14.530 . we used the background equations ( 12.19^ in 
the slow-roll approximation of the subsection 12. 2[ The conformal time r has the following relation 



dt 
a 



^H 



aH / a?H 



-da 



aH 



1 / a'' Ml 



6a2/72 



where the relation 



.H 



H ^ 



«2M2 



(4.55) 



(4.56) 



was used. Here we keep the leading order in the slow-roll parameters, e and rj, and a small gradient 
constant a. We will discuss the range of a in our approximation later. Then, from (14.551) . we 
obtain 



^2 ^ l+2£ ^ «2M| 



(4.57) 



r- 3 ' 

where "H = a' /a. 

Assuming that the slow-roll parameters are constants, we have the exact solutions for the 
equations of (I4.53p . 



V, 



Uk 



'— r 



C,kH^^){-k,T)+d,kH^^^{-hT) 
C2kH^^X-hT)+d2kHJ^J{-hT) 



ei{k), 



e2{k), 



(4.5^ 



where HJ^ (z) {i = 1,2) are the first and second kinds of the Hankel functions, and ei(/c)'s are 
independent Gaussian random variables satisfying 



(e,(fc)) = 0, {e,{k)e,{k')) = 5,,5{k-k'). 
Here, the angled brackets denote ensemble averages. 



(4.59) 
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In order to determine the integration constants Cik and dik-, we have to fix the initial vacuum 
state at /cr — )■ — oo. As we see in ( 14.53^ and (14.58^ . the two differential equations and the 
corresponding solutions have the same forms. Due to the overall normalization of u mode given 
in (13 .29 p . Vk and Uk modes have the same normalizations. Therefore, we can concentrate on the 
solution of Vk only and apply the result to the case of Uk- 

The mode for Vk = Vkei{k) satisfies the normalization condition, Vkvl — v^Vj^ = i which is 
originated from the quantization condition of the mode. Due to this normalization, we have the 
relation between the integration constants, 

\cik\^ - \dik\^ = ^. (4.60) 

We adopt the Bunch-Davies vacuum for the initial perturbation mode at r ^ by taking the 
positive energy mode. Then the initial mode is given by 

Vk{T) = ^e-^'^^e,{k). (4.61) 



This mode corresponds to the choice of coefficients. 



ci, = e^('^+^)^, d,k = 0. (4.62) 

Using the relations fl4.6ip and fl4.62p . and adopting to the case of Uk as well, we have the exact 
solutions for Vk and Uk, 



Vkir) = ^e^(-^+^)v/^/7«(-A:,r)e,(fc) 



Wfc(r) = ^e'^(^^+-^)^/^H^^J{-hT) e^ik). (4.63) 

As we see in the initial mode fl4.6ip . the effective wave number ki should be real to have a 
well-defined quantum state. This means that the square of the effective wave number kf should 
be non-negative. Therefore, the wave number k should be constrained from fl4.54p as 

k' > kl, ^ —^. (4.64) 

As we see in f l4.53p . the initial quantum fluctuation modes of Vk with |A;| < |A;min| are exponentially 
growing or decreasing and so they break the unitary symmetry of the quantum system in the deep 
inside of the horizon. That is, these modes are not allowed as initial quantum modes. Therefore, 
there is an observational lower bound of the comoving wavenumber in our model. In the large 
scale limit A;i.2|t| <^ 1, the modes in (14.631) can be approximated as 



.iu„,_3r(/ii) 1 



1 



Vk{r) - e^l'^^+5^2^^-2 -^_=(-A;ir)2-^i e,{k), 

-'- 12'' V 2fci 



l^_„_ 3 i /i9 i , , ,1 



.3r(/i2) 1 

.2' 



Ukir) ^ e^l'^^+5j2^'^"2 -^ (-A;2r)5-^2 e2{k). (4.65) 
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In the slow-roll approximation with .^ = and a small nonvanishing value of a, the leading 
contributions for the curvature and isocurvature perturbations in (13.43^ are given by 

From these expressions, we obtain power spectrums for TZ and 5, 



^^ ,^^.v fc^a^ 



P5(fc) = 77^(55*) 



(^^ 



^(^g:) + ^ ((Q.g;) + (Q.g;)) 



(4.68) 



27r2 ^ ' 27r2e2 

Using the relations of (14.521) and (I4.65p . we obtain the power spectrums at the horizon crossing 
point {k = aH) as 



^-- - (l8^ j ?II (1 + 2C^ - 32. - 18., - ^^ j . (4.69) 

where the subscripted asterisk indicates the value at the horizon crossing point and C = 2 — In 2— 7 
with the Euler-Mascheroni constant 7 ^ 0.5772. In the computations of power spectrums in (14. 69 p . 
we have used the relations, 

(Q.Q;)* = ^(V.V*). ^ g (^1 + (6C - 8)6 + (2 - 2C)v - ^ 



(QuQl) = ^{UM)* ^ ^ (^1 + 2(C - 2)e 



The expansions of power spectrums in (I4.69P are only valid in the range of a, 

Mia" < e^e. (4.71) 

In a — )■ limit, the power spectrums Vtz and Vs at the horizon crossing point in (I4.69P become 
those of the single field case. Finally, the spectral indices for TZ and S in the large scale limit 
(^i,2|t| <^ 1) are given by 

rfln^T^ _^ ^ , 2a^M^ 



nTi-l = —77-^- ~ 2?7 - 6e , ^^^ 

^s-l^ -n-7^ ^ -2 + 4ry- 10e+ — — ^. (4.72) 



dink ek"^ 
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If we take (14.711) into account, we find the last terms in n-ji — 1 and n^ — 1 are of the order of e 
and e^, respectively. These imply that while the curvature perturbation is nearly scale invariant, 
the isocurvature perturbation is proportional to the inverse square of the comoving wavenumber. 
The running of the scalar spectral index is given by 



d\nk ek^ 



+ 0{e\er^). (4.73) 



The leading term becomes of the order of e from (14.691) and (14.711) . Recent Planck [3] and WMAP 
nine-year |1] data show 

rin = 0.9603, -^ = -0.0134 (Planck), (4.74) 

dmk 

nn = 0.9608, -^ = -0.019 (WMAP - 9yr) (4.75) 

with the pivot scale ko = 0.05 Mpc"^ (Planck) and ko = 0.002 Mpc"^ (WMAP), respectively. 
Since a^ > and e > 0, the sign of the running of the spectral index in (I4.73P is consistent with 
the observations. However, to give some restriction on the value of a we need more investigations 
using a specific potential of the single scalar field. We leave these issues for future work. 

4.2 Tensor mode 

Using the relations (14.561) and (I4.57p . we obtain 

Then the linearized perturbed equation (I3.47P for the tensor mode, up to the leading order in the 
slow-roll parameters and the gradient constant a, can be written as 



u'' 



k 



11^ - 1 

2 _ PT 4 



ut = 0, (4.77) 



where 



4 = ^ 6"^ ' /^^ = 2 + '• ^^-^^^ 

Interestingly, the tensor mode has the same lower bound of the comoving wavenumber with that 
of the scalar mode in (14.641) . 
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Choosing the initial condition with a positive frequency mode, just hke fl4.6ip in the initial 
state of the scalar mode, we obtain the following solutions in the large scale limit (A;t|t| <^ 1), 



where A{e) = 1 — eln(— fcr) + e(2 — 7 — ln2) and cxik) is an independent Gaussian random 
variable and has the same properties with those defined in fl4.59p . 

With these solutions, the power spectrum for the tensor modes hij is obtained as 

where we used h\ = -^;jj-ut and the factor "2" comes from the two polarization states. We also 
used the quantities at the horizon crossing point {k = aH), 

12A;2 



;-A;Tr)"|,^ l + e + 



{UTU*^). ^ ^ (1 + (2C - 4)e - ^) . (4.80) 

The spectral index for tensor modes and tensor-to-scalar ratio are given by 

r = :^ ~ 16e. 1 + (4 - ^C)e + (2C - 2)r/ + "^—^ . (4.82) 

Vn \ k^e ) 

5 Conclusion 

In this paper, we studied the cosmological properties of a multi-field inflation model in which we 
considered the single field and an additional triad of scalar fields with a non-minimal coupling 
which is dependent on the single scalar field. 

As a special background solution, we considered spatially dependent linear profiles for the triad 
of scalar fields, while the single scalar field depends on the cosmic time only. Due to this linear 
property, all components of the energy momentum tensor evolve homogeneously and isotropically 
on the background FRW metric. For this reason, we can regard the background evolution as that 
of the single field inflation in the presence of a matter with an equation of state w = —1/3, which 
spreads to the whole space. However, there appeared one more scalar degree of freedom from the 
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triad of scalar fields in the perturbation level. Therefore, our model is different form the single 
field inflation model in some cosmological background. 

For the exponential-type potential, we found an exact solution describing the power-law infla- 
tion, as in the case of the single field inflation model. We also investigated the properties of the 
background evolution under the assumption of the slow-roll inflation. The background behaviors 
of the slow-roll inflation are determined by the two parameters ^ and a which are associated with 
the non-minimal coupling and the gradient constant of the triad of scalar fields, respectively. We 
have drawn contour lines of the e-foldings for the quadratic potential on the parameter space of ^ 
and a in Fig.l. 

We also investigated the properties of the scalar and tensor perturbations for the slow-roll 
inflation with vanishing ^ and small a. In the scalar fluctuations, there are two physical modes, 
Vfc and lAk originating from the scalar field and the triad of scalar fields respectively in the spatially 
flat gauge. The two scalar modes are decoupled from each other and have the form of Sasaki- 
Mukhanov equation in the leading contribution of a. We found that there exists an observational 
lower bound of the comoving wave number, |/cmin| ~ aMp, to have a well-defined initial quantum 
state for Vk with a positive energy eigenvalue. We also obtained the power spectrums and the 
corresponding spectral indices for the comoving curvature and isocurvature perturbations. If we 
compare our results for the power spectrum with those of the single field inflation model, the 
a-dependent terms are always proportional to a'^/k^. Due to this reason, the spectral indices for 
the curvature and isocurvature perturbations have the contribution of a^ /k"^ terms as well. As is 
already well-known, the isocurvature perturbation in the single field inflation model is vanishing in 
the large scale. However, the power spectrum for the isocurvature perturbation is non-vanishing 
and has contributions from both V^ and Uk modes in our model. We also implemented the same 
analysis for the tensor perturbation. We found the same lower bound of the comoving wavenumber 
in the tensor mode. 

In the case of non-vanishing ^, Vk and Uk modes are coupled with each other. Then it is 
difficult to obtain analytic solutions, unlike the case of ^ = 0. As we saw in the analysis of 
the background evolution, the value of ^ changes the behaviors of the cosmological evolution 
significantly. Therefore, we can expect that the value ^ may have some important roles in the 
perturbations as well. We leave the analysis of the nonvanishing ^ for future work. 

One important difference between the single field model and ours is the existence of the isocur- 
vature perturbation in our model, which originates from the spatial gradient of the background 
solution. One may also test whether the nonvanishing spatial gradient can be a source of the 
cosmological non-Gaussianity and compare the result with the recent observational data. 
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